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Abstract 

A beautiful understanding of the smallness of the neutrino masses may be ob- 
tained via the seesaw mechanism, whereby one takes advantage of the key qualitative 
distinction between the neutrinos and the other fermions: right-handed neutrinos 
are gauge singlets, and may therefore have large Majorana masses. The standard 
seesaw mechanism, however, does not address the apparent lack of hierarchy in the 
neutrino masses compared to the quarks and charged leptons, nor the large leptonic 
mixing angles compared to the small angles of the CKM matrix. In this paper, 
we will show that the singlet nature of the right-handed neutrinos may be taken 
advantage of in one further way in order to solve these remaining problems: Unlike 
particles with gauge interactions, whose numbers are constrained by anomaly can- 
cellation, the number of gauge singlet particles is essentially undetermined. If large 
numbers of gauge singlet fermions are present at high energies - as is suggested, for 
example, by various string constructions - then the effective low energy neutrino 
mass matrix may be determined as a sum over many distinct Yukawa couplings, with 
the largest ones being the most important. This can reduce hierarchy, and lead to 
large mixing angles. Assuming a statistical distribution of fundamental parameters, 
we will show that this scenario leads to a good fit to low energy phenomenology, 
with only a few qualitative assumptions guided by the known quark and lepton 
masses. The scenario leads to predictions of a normal hierarchy for the neutrino 
masses, and a value for the |m ee | mass matrix element of about 1 — 6 meV. 



1 Introduction 



The quark and lepton masses and mixing angles provide one of the fundamental mysteries 
of the standard model of particle physics. They hint at a potentially deep underlying 
structure, while also seeming sufficiently random so as to defy straightforward explanation. 
There are a number of qualitative features of these parameters which one would like to 
understand!]] 

1. Why are the neutrino masses roughly ten orders of magnitude smaller than the 
masses of the quarks and charged leptons? 

2. Why do the quark and charged lepton masses have significant hierarchies of about 
five orders of magnitude? 

3. Why is the CKM matrix approximately equal to the identity when the up and down 
quarks are both ordered by mass? 

4. Why are the mixing angles in the lepton sector fundamentally different from those 
in the quark sector, with two angles close to maximal? 

Of these questions, it is the first which has lent itself most easily to explanation. The 

key observation is that right-handed neutrinos, unlike any of the other standard model 

fermions, have no known gauge interactions. As a result, no symmetry forbids them from 

obtaining Majorana masses, which may be many orders of magnitude larger than the 

weak scale. After integrating out the right-handed neutrinos from the theory, one obtains 

effective operators of the form (HL) 2 /Mr, where H is the Higgs field, L is a lepton doublet, 

and Mr is the right-handed neutrino mass scale. This yields light Majorana neutrinos with 

masses of order v 2 /Mr, with v = 174 GeV being the Higgs vacuum expectation value, and 

with a possible additional suppression from Yukawa couplings pQ. Appropriately small 

masses are then obtained for M of order 10 14 — 10 16 GeV, perhaps related to the scale 

of grand unification. This picture is so simple that it is now taken almost for granted, 

although it should be kept in mind that experimental confirmation is still lacking. 

1 The strong CP problem may or may not have a solution related to the other flavor mysteries, and 
we will not concern ourselves with it in this paper. We may rely, for example, on a separate solution 
involving an invisible axion with a large Peccei-Quinn symmetry breaking scale. 
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The remaining questions on our list, on the other hand, have no such obvious inter- 
pretations, although a large variety of proposals have certainly been put forth. In many 
models, one assumes that the quarks and leptons are charged under a variety of possi- 
ble flavor symmetries, with the Yukawa couplings proportional to appropriate powers of 
symmetry breaking spurions. While such models can produce successful phenomenology, 
the choices for the flavor symmetry transformations, as well as the associated spurion 
structures, can seem ad hoc, and it is perhaps fair to say that no particular model stands 
out as being especially compelling. 

Another possibility, which runs somewhat counter to the flavor symmetry perspective, 
is that the various Yukawa couplings are determined in some more complicated way by the 
fundamental theory, with their low energy values not representing any essential pattern. In 
this scenario, one simply parametrizes the Yukawa couplings via a hypothesized statistical 
distribution, and checks the phenomenological consequences. Such distributions may 
arise, for example, via the landscape of string theory, or may perhaps be manifestations 
of our own ignorance of the fundamental theory. At the very least, the extent to which 
such an approach is viable is an important issue for determining the merit of the flavor 
symmetry paradigm. Fortunately, the standard model contains a large number of Yukawa 
couplings, and these may then serve as a reasonable statistical sample from which to draw 
proposals for the Yukawa distributions. A key issue with this approach, however, comes 
from question number (4), above. Naively, it seems difficult to allow for a single form 
of probability distribution to describe all of the standard model Yukawa couplings. The 
quarks and charged leptons seem to prefer a probability distribution which scans roughly 
evenly over many orders of magnitude in order to yield large mass hierarchies and the 
small CKM mixing angles [2]. The neutrinos, on the other hand, would seem to prefer a 
more degenerate - or "anarchical" type of scan in order to obtain the required large mixing 
angles [31 SI E] - While we currently have no way to derive these probability distributions 
from fundamental theory, and it is certainly possible that the neutrino Yukawas are simply 
determined in a special way, different from the other fermions, a unified description would 
seem more attractive, and would perhaps lend more weight to this approach. 

In this paper we will show that the seesaw mechanism alone may be what is distin- 
guishing the neutrinos from the other fermions, both in terms of their masses, and also in 
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terms of their large mixing angles. In particular, we will show that it is indeed possible to 
make use of the same general form of distribution for all of the Yukawa couplings of the 
standard model, and still produce appropriate phenomenology. We will consider scans 
with large hierarchies, which seem naively most suited to the quarks and charged leptons. 
As in the standard seesaw, for the neutrinos we shall take advantage of the fundamental 
physical feature which distinguishes them from the other fermions; the fact that the right- 
handed neutrinos are gauge singlets. As a result, not only do the right-handed neutrinos 
have the possibility for large Majorana masses, but moreover, due to the absence of any 
anomaly cancellation requirements, we do not know how many of them there are. In- 
deed, constructions in string theory often produce large numbers of singlet fermions with 
Majorana masses close to the GUT scale after the compactification of extra dimensions 



In the presence of a large number of singlets, the low energy neutrino mass matrix will 
be realized as a sum over many distinct hierarchical numbers with the largest Yukawa 
couplings dominating. This generically results in a washed out hierarchy, and potentially 
large mixing angles. It follows that a random hierarchical scan of Yukawa couplings is in 
fact a good fit for the observed masses and mixing angles of the standard model, so long 
as the possibility for a large number of right-handed neutrinos is taken into account. 

The reason that many singlet states often arise with masses close to the GUT scale 
in string constructions is actually reasonably generic: After compactifying d extra dimen- 
sions in a string model, the resulting Kaluza Klein mass scale is related to the effective 
four-dimensional Planck scale Mp\ through the relation 



where M s is the string scale. In order for the compactification geometry to be reasonably 
described by a classical gravity picture, it is required that Mkk be at least parametri- 
cally smaller than M s . From the above relation, and considering the 6 available extra 
dimensions of string theory, we thus see that it may be reasonable to expect the Kaluza 
Klein scale to be roughly a few orders of magnitude smaller than MpiJ^] Note that it is 
2 For more general field theory discussions, see also [51 [TU], 

3 It is not necessarily required, however, that the compactification of the extra dimensions have a well 
defined classical gravity description. This is perhaps the main caveat to the present argument. 




(i) 
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also possible that grand unified symmetry is broken via a mechanism directly related to 
the compactified extra dimensions, further motivating the presence of the KK scale at 
close to 10 14 — 10 16 GeV. In any case, the point then is that the KK scale is a natural 
scale for the appearance of a large number of gauge singlets which may then serve as 
right-handed neutrinos; these singlets might be born out of KK towers themselves, as was 
the case considered in [7], or they might be associated with moduli fields, stabilized by 
masses close to Mkk due to fluxes in the compactified dimensions |8] Q For our purposes, 
the main point is that the number of relevant singlets with masses close to the needed 
seesaw scale may easily number in the tens or hundreds. Note that we will not attempt to 
construct an explicit top-down model for Yukawa couplings in this paper, but will simply 
take a phenomenological point of view based on the known properties of the standard 
model masses and mixing angles. We will leave top-down model building as a possible 
subject for future work. 

The outline of our paper is as follows: In section 2 we will describe our framework in 
more detail, discussing the various qualitative features we require, including the properties 
of the high energy Majorana mass matrix. In section 3, we will give an example of a 
particular universal Yukawa coupling distribution which gives a good phenomenological 
fit for the charged fermion masses and mixings, and show how our mechanism then also 
leads to good phenomenology in the neutrino sector. Additionally, we will show that our 
setup leads to the predictions of a normal hierarchy for the neutrino masses, and a value 
of |m ee | - the neutrino mass parameter relevant for neutrinoless double beta decay - of 
between about 1 and 6 meV. We will summarize and discuss future directions in section 
4. 

2 Framework 

The reason that many right-handed neutrinos, which we shall denote generically as u' R s, 
can wash out hierarchy in the low energy neutrino mass matrix m is simple. This matrix 

4 In this latter case note that the singlet masses are actually given by ~ M^ K /Mg, with a further 
small suppression below the KK scale. 
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takes the form 

mj^v^YikM^Yfl, (2) 

Ik 

where, with N right-handed neutrinos, Y is the 3 x N Yukawa matrix, M is the N x N 
Majorana mass matrix, and v is the Higgs vacuum expectation value. Taking the typical 
largest Yukawa coupling size to be of order one, for example, we can see that, with 
sufficiently many z/r's, each matrix element in m will generically obtain some number of 
large contributions. Even if the original Yukawas were hierarchical, the m^'s then come 
out to be roughly of the same order, with the hierarchy having been lost (note that, 
assuming arbitrary signs or phases for each term in the sum, one obtains an enhancement 
to the overall scale of m by a factor which scales as \N). 

The mechanism we are discussing here is fairly general, and a specific structure or 
distribution for the Yukawa couplings and Majorana masses is not necessary. We do, 
however require the following conditions to be satisfied: 

• There must be a reasonably well defined upper-bound to the neutrino Yukawa cou- 
plings. Such an upper-bound could be set by perturb at ivity or by some other fun- 
damental physics, and should be reasonably independent of which neutrino field a 
given Yukawa coupling is associated with. 

• Enough right-handed neutrinos must be present so that each left-handed neutrino is 
expected to have ^ 1 Yukawa couplings within a factor of a few of the upper bound. 
It can be seen by inspection of equation [2] that this would result in each low energy 
neutrino mass matrix element expecting to receive at least one "large" contribution. 
If the Majorana mass matrix has hierarchical eigenvalues, then this condition applies 
to the lightest set of the z/_r's, which will give the dominant contribution to the seesaw. 
Amongst the standard model quarks and charged leptons, we have one order one 
Yukawa coupling amongst 27 (associated with the top quark), and so with this as 
a guide, we anticipate that we may require roughly ^ 30 right handed-neutrinos for 
our mechanism to work. 

• We require a Majorana mass matrix which mixes together the right-handed neutrino 
fields. One simple way this might be realized is if the physics determining the 
Majorana masses is distinct from that determining the Yukawa couplings (in which 
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case the unitary matrix which diagonalizes the Majorana masses can be taken to be 
distributed according to the Haar measure, as will be discussed in section 3), but 
essentially any sufficiently non-diagonal structure should be viable. 

The need for this last condition can be seen as follows: Suppose that the Majorana mass 
matrix were diagonal. In that case there would be a tendency for the diagonal elements 
of the low energy mass matrix to be larger than the off-diagonal elements, resulting in 
small mixing angles. Let us compare for example the 1, 1 element to the 1, 2 element: 

"n~Ett ( 3 ) 

k 

m la ~X) M *fc ly ifc y »- ( 4 ) 

k 

For any given Yukawa coupling y, there is in general a much higher chance that y 2 will 
be near its upper bound, than the chance that two separate couplings y\ and yi will both 
be large and lead to a large product y\yi- It is thus clear from equations [3] and [4] that a 
diagonal Majorana mass matrix will tend to give poor phenomenology in our scenario]^] 

It should be clear that we certainly do not require that all of the Yukawa couplings 
of the quarks and leptons share a universal distribution. We will, however, concentrate 
on such distributions in this work, since in this way we may check that it is actually the 
many z/r's which are resulting in large mixing angles, rather than any fundamental lack 
of hierarchy assumed for the neutrino Yukawa couplings. Note that for this reason, there 
is no intrinsic obstacle to realizing our scenario in a supersymmetric context (in which 
tan/3 would modify the hypothetical quark and charged lepton distributions somewhat), 
or even in the context of a grand unified theory. 

Note that a wave-function overlap picture (see e.g. [3 El EH 112]) for obtaining hierarchy 

in the Yukawa couplings would not work very well for our scenario, since in such cases 

the couplings for a given field tend to all be correlated in size. This is the key manner in 

which our scenario with many right-handed neutrinos differs from those already appearing 

in the literature. Constructions with compactified extra-dimensions may still be used 

5 Note that, even with a Majorana mass matrix with substantial off-diagonal components, squared 
Yukawa couplings will still only contribute to the diagonal elements of my. With many right-handed 
neutrinos, however, the large number of terms (of order N 2 ) which are products of separate Yukawa 
couplings actually give the dominant contribution. 
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to motivate the presence of the many ur states near the GUT scale, but we require, 
for example, extra-dimensional wave function profiles to be fairly flat, with the Yukawa 
hierarchy generated by fundamental physics in some other manner. 

In the next section we will show our mechanism at work quantitatively in a specific 
example: We will take a Yukawa coupling distribution fit to the known properties of 
the quarks and charged leptons, and work out the phenomenological consequences as a 
function of the number of VrS. We will also show that our scenario (as well as any 
scenario with a roughly anarchical low energy neutrino mass matrix) leads to a prediction 
of |m ee | ~ 1 — 6 meV. 

3 An Example 

Fitting to the charged fermion masses and CKM angles 

Without knowledge of the high-scale theory it is of course impossible to know the fun- 
damental distribution of the Yukawa couplings. However, existing data can serve as a 
guide; the best fitting probability distribution for the charged fermions has been exam- 
ined in some detail in |2], which we use as a starting point for our study. For power 
law distributions with minimum and maximum cutoffs, they found that the quark and 
charged lepton masses are best fit by a distribution that is very close to scale invariant, 
p(m) ~ 1/m. However, we expect that the Yukawa couplings, not the masses, are the 
fundamental parameters of the underlying theory. Because of mixing effects, the form of 
the Yukawa probability distribution is not directly transferred to the masses. Considering 
distributions of the form 

p{y) ~ ^ (5) 

at the GUT scale within the standard model, the authors of |2] found that values of 5 
just below 1.2 give the optimal mass distribution, with 1.2 x 10 -6 < y < 1. 

The Yukawa couplings also determine fermion mixing, so we wish to ensure that our 
choice of distribution also leads to reasonable values for the CKM elements. In fact, 
the hierarchical structure considered here naturally leads to small mixing angles. Note 
however, that there is a caveat; without any flavor symmetry to correlate the Yukawa 
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couplings of a given field, we must rely on a minor coincidence to ensure that it is actually 
the mass ordered basis in which the CKM matrix appears approximately diagonal. The 
probability that this could have occurred by chance is 1/3! ~ 17%, and we will take this 
as a starting assumption, considering only CKM matrices in which the largest element of 
each row and column lies on the diagonal in the mass ordered basis. For the distribution 
considered in [2], with p(y) ~ 1/y 116 , while it was found that reasonable masses and small 
CKM mixing were obtained, in some cases the observed CKM off-diagonal elements were 
larger than the mean predicted values, lying on tails of their distributions. This can be 
seen from the plot of the \V US \ distribution in Figure [TJ 
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Figure 1: |V^ S | for 5 = 1.1 for matrices with mass eigenvalues within a factor of VlO of 
GUT scale values. The vertical red line represents the experimental value of \V US \. The 
maximum possible value is l/\/2 due to our requirement that the diagonal entries be the 
largest. Inset: Cumulative probability distribution function for 8 = 1.0 — 1.3. Vertical 
dotted lines indicate the median values of the distributions, and horizontal dotted lines 
indicate the fraction of matrices with \V US \ less than the experimental value. 



To ensure that we consider a distribution which gives reasonable quark mixings, we 
have generated a large number of up and down quark mass matrices, by scanning according 
to equation [5] for various values of 5 and assigning each element a random phase. We 
then consider those which give a mass eigenvalue for each quark which lies within an 
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Table V. Median values of CKM element magnitudes for different values of 5 along with 
their measured values |15j . The rightmost columns indicate the percent of matrices which 
give a given CKM element smaller in magnitude than the measured value. Restrictions 
on the mass eigenvalues are as in Figure [TJ 

order of magnitude of the measured value, i.e. in the range (m q /y/lO, m q y/l0)f\ Finally, 
we diagonalize the mass matrices with eigenvalues ordered by mass and determine the 
corresponding CKM matrices, eliminating any with an incorrect generation structure, 
as discussed above. As we see in Figure [TJ there is a sharp peak at small mixing, so 
to get a measure of how consistent each scenario is with the data, we characterize each 
distribution by its median as well as the probability of finding a value smaller than the 
experimentally observed value. We will use the 68% and 95% levels of the latter to define 
ler and 2a ranges, as in [2]. This is illustrated for \V US \ in Figure [TJ (inset), which shows 
some preference for smaller values of 5 - smaller 8 leads to larger mixing angles. To get an 
overall measure of the CKM matrix, we consider the three off-diagonal elements shown in 
Table [TJ The measured values of \V US \, \V ub \, and \ V cb \, all show improvement of agreement 
as 5 is decreased, and fall within our la criterion around 5 = 1.1. We use this, along 
with the mass only best fit of 5 ~ 1.2 ± .1 from [2], in our choice of distribution. While 
one could perform a detailed statistical analysis to determine the best fit to all measured 
parameters, the 5 = 1.1 case gives rise to both reasonable charged fermion masses as well 
as CKM angles, so we will take it as our example distribution for studying the neutrino 
sector. 

6 Here and throughout this paper we use charged fermion mass eigenvalues evaluated at the GUT scale 
in the standard model [T3] , though [2\ demonstrated that the best fit for the distribution feels little effect 
from the running. The experimental values for the CKM elements we cite do not reflect any running, 
but this effect is small, no more than ~ 15%, and is negligible for V us fH] . 
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The Majorana distribution 

For the high energy Majorana mass matrix, for which we have no guidance from data, we 
take a simple ansatz in accord with the discussion of section 2; we suppose that the physics 
determining the Majorana mass matrix is independent from that determining the Yukawa 
couplings. This assumption then implies that the unitary matrix U which diagonalizes 
M = UDU T may be taken as distributed according to the Haar measure Haar 
distributed unitary matrices are easily generated via previously established algorithms 
[To] . Note that we must further impose a distribution on the mass eigenvalues appearing 
in the diagonal matrix D; for simplicity we choose the case of a linear distribution between 
zero and Mr, the Majorana mass scale, which we fix by requiring that we match the value 
for the larger measured neutrino mass difference, Am^. 

Probability for large PMNS mixing angles 

Given our chosen statistical distributions, we next generate a large number of PMNS 
matrices for various numbers of right-handed neutrino species between N = 2 and N = 
100. We first generate a3xiV Dirac mass matrix Md = vY, where the entries of the 
Yukawa matrix Y scan between 1.2 x 10~ 6 and 1 with probability density p(y) ~ y~ 1A and 
have random phases. We then generate a Haar distributed N x N Majorana mass matrix 
with linearly scanning eigenvalues, to form a neutrino mass matrix m = MdM~ 1 M'£ = 
U^D^Uj, where D u is a diagonal matrix containing the neutrino mass eigenvalues, and 
U v is unitary. We generate 3x3 charged lepton matrices Mi with the same distribution 
as the Dirac matrices for neutrinos. To more appropriately model the true charged lepton 
mass matrix, we only take matrices which agree to within a factor of \/l0 of the GUT 
scale charged lepton masses, as described previously for the quarks. We diagonalize 
Mi = UiDiV?, from which we can generate the PMNS matrix Upmns — UjU u . This 
can be parameterized, with unphysical phases rotated away, in the standard fashion with 

7 The Haar measure /iona group G satisfies n(gE) = n(E) for every g £ G, and is uniquely denned 
for G=U(N). 
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three mixing angles (6^) and three CP phases (5, 021, «3i' 



U. 



PMNS 



C12C13 

— S 12C23 — Ci2S23Si3e* 
iS 

Cl2C23 s 13 e 

21 „- "31 



S12S23 

x diag(l,e*^ i ,e i ^ i ) 



S12C13 

iS 

C12C23 — Sl2S23 s 13 e 
— C12S23 — Sl2C23 s 13 e * 




(6) 
(7) 



where s^- = sin(^) and qj = sin(6 l ^)|^] Experimentally, there are two large mixing angles, 



0.312ig;g^ and maximal 
mixing for sin 2 (#23).B.F = 0.521q'q7 (lcr)[17j. The third mixing angle is known to be small, 
with a best fit of sin 2 (^ 13 ) Bi r = 0.013+q;qq5 at la [T7j, and there is growing evidence for a 
nonzero mixing jT8j [19] . 
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Figure 2: Fraction of matrices with sin 2 (26 l )i ar g CSt > 0.98 and sin 2 (2#) 



/next largest 



> 0.84. 



To see how common such large mixing angles are in our scenario, for various num- 
bers of right-handed neutrinos, we determine the fraction of matrices having both one 

8 Note that in this paper we ignore running effects due to the neutrino Yukawa couplings. Even at 
and above the seesaw scale these are expected to be small for the distributions we consider, since we 
take all of our Yukawa couplings to be smaller than 1, with only a handful approaching this upper 
bound. If the bound on the Yukawas were taken larger, then renormalization group running could have 
important, albeit model dependent effects. For example, with low energy supersymmetry, depending on 
the superpartner spectrum, it is possible for running at the seesaw scale to contribute to slepton mass 
splittings and flavor changing neutral currents [9]. If only the standard model were assumed valid up 
to energies above the seesaw scale, then it is possible to obtain large radiative corrections to the Higgs 
quartic coupling and thereby large, experimentally excluded, values for the Higgs mass. Note that this 
latter conclusion would be avoided in the presence of high scale supersymmetry. 
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near maximal mixing angle and one large mixing angle. Because maximal mixing, with 
sin 2 (2$) = 1, is a special point, we look for cases which have at least as much mixing 
as the la experimental bounds, requiring that one angle satisfies sin 2 {26) > 0.98 and 
another satisfies sin 2 (2#) > 0.84. The results are shown in Figure [2J from which we see 
a clear indication that as the number of right-handed neutrinos increases, so too does 
the likelihood of obtaining large mixing angles - as expected for the reasons laid out in 
Section 2. This effect is further illustrated in Figure |3j where we see the shift to larger 
mixing angles as N increases. 
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sin 2 (26>) 

Figure 3: Distribution of mixing angles. The three different bands represent the largest, 
middle, and smallest sin 2 (29). 



Other parameters 

While the absolute masses of the neutrinos are not well measured, oscillation experiments 
give us a good measure of their mass squared differences, with a best fit of Am^ = 
7.59±o; 2 g x 10~ 5 eV 2 and Am 2 x = 2.50to;?g x 10~ 3 e ^ 2 (assuming a normal hierarchy, with 
comparable values for an inverted hierarchy) [T7j. To see if our construction accommodates 
this small but non-trivial hierarchy, and to determine whether there is a preference for a 
normal or inverted structure, in Figure [4] we consider the ratio of neutrino mass squared 
differences, which we plot as log 10 Amf^/Am^. Here we label the masses such that 
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m3 > m 2 > mi, so that this quantity is positive for a normal hierarchy and negative 
for an inverted onej^] Observed masses give a value of about ±1.5. We see that for 
large N, the masses are much less hierarchical, and easily accommodate the observed 
values. Furthermore, we see an overwhelming preference for the normal hierarchy, which 
in particular justifies our use of the associated mass and mixing angle measurements in 
later parts of this section, 
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Figure 4: Ratio of mass squared differences log 10 Am^/Am^ for N = 3, 10, 30 and 
100. Here we choose the convention m 3 > m 2 > mi, so that positive (negative) values 
correspond to a normal (inverted) hierarchy. 



Having seen that the mixing angles and mass splittings observed in nature are in- 
creasingly typical as N increases, we wish to look at other properties of viable matrices 
produced within our framework. To select cases close to reality, we consider only matrices 
which satisfy: 0.28 < sin 2 (# 12 ) < 0.35; 0.41 < sin 2 (# 23 ) < 0.61; 29.1 < Am^/Am^ < 
35.6; and 0.004 < sin 2 (#13) < 0.028, which come from best fit 2a bounds [T7j. In Figure 
[5| we show the distribution of sin(6 , 13 ), subject to the large angle and mass constraints, 
and find that there is some tension with the best fit, which at 2a corresponds to about 
9 Note that for an inverted hierarchy, our labeling is non-standard. 

10 The reason our scenario strongly prefers a normal versus an inverted hierarchy is that the reasonably 
large observed ratio of solar and atmospheric mass squared differences necessitates that either the heaviest 
(normal hierarchy) or the lightest (inverted hierarchy) of the neutrinos is a mild outlier. Having the 
heaviest neutrino as the outlier in our scenario is much more probable, since this requires fewer outlying 
elements in our typically degenerate mass matrix. 
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0.06 < sin(#i3) < 0.17. On the other hand, after fitting successfully the other parameters 
of the neutrino mass matrix, the probability of obtaining one mild outlier is not too small. 
This can be seen on the right side of Figure |5j which gives the fraction of matrices at or 
below a given value of sin(#i3) - there is significant variation across the experimentally 
preferred region (in grey), reaching around 10% near the upper 2a bound. Our scenario 
clearly prefers a nonzero value for sin(#i3), for which there is growing evidence, and a 
significant range of allowed values would be unexceptional. The situation for #13 in our 
scenario is ultimately similar to that of neutrino anarchy, for which a global statistical 
analysis was performed and found good agreement with the data [5] . Performing a similar 
analysis for the present case is beyond the scope of this work, but we would expect to 
find similar results. 



2.0 rx 




sin (613) sin(<?i 3 ) 

Figure 5: Left: sin(# 13 ) distributions among matrices satisfying 2a large angle and mass 
constraints described in the text, for N = 30 (blue, solid) and iV = 100 (red, dashed). 
Right: The corresponding cumulative probability distribution functions. The grey shaded 
area indicates the 2a best fit region. 

We next consider the Majorana mass scale, which we determine by requiring that the 
mass splitting Am^ matches its experimental value. We plot the resulting distribution in 
Figure [6j where we see that the distribution rises and sharpens at large N, with typical 
values around 10 15 — 10 16 GeV, suggestive of the framework discussed in Section 2. 

Neutrinoless double beta decay provides us with an experimental probe of Majorana 
interactions of neutrinos, and the figure of merit for such experiments is given by m ee , the 
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Figure 6: Majorana mass scale distribution among matrices satisfying 2a angle and mass 
constraints described in the text, for N = 30 (blue, solid) and N = 100 (red, dashed). 

upper left entry of the neutrino mass matrix in the diagonal charged lepton basis, 

m ee = U^rrii, (8) 

i 

where U is the PMNS matrix and rrii are the masses of the light neutrinos. We show 

the prediction for m ee in Figure [7j again using only matrices lying sufficiently close to 

observation, and with the overall scale set by the true value of Am^. The 68% central 

regions of the distributions are .0015 - .0049 eV for iV = 30 and .0018 - .0062 eV for 

iV = 100, which seems typical for theories in which the low energy neutrino mass matrix is 

composed of relatively degenerate numbers - for comparison, we also show the distribution 

we get by applying the same procedure to matrices with random 0(1) parameters, i.e. 

anarchical matrices]^] Unfortunately, the predicted values for m ee are beyond the reach 

of current experiments. On a 10 year or longer timescale, it is possible that neutrinoless 

double beta decay searches, such as a future iteration of EXO [20J, or observations of 

large scale structure j2TJ E2] may reach the required sensitivity. On shorter timescales, 

n Here we generate anarchical matrices with basis-independent scanning following the prescription of [3]: 
We generate a complex 3x3 Dirac matrix, Mp, and complex symmetric 3x3 Majorana matrix, M, each 
with entries scanning between and 1 with random phases, subject to the basis-independent boundaries 
Tr(MU/), Tr(MjjM£>) < 1. We then form our light neutrino mass matrix, m = AM D M~ 1 M'£, where A 
is chosen to obtain the correct value of Am^. Because of the basis independence, we are free to choose 
the charged lepton matrix to be diagonal. 
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it is interesting to note that our framework could be falsified if neutrinoless double beta 
decay were to be observed with a large measured value of m ee . 

f 300 1 : 

> 1 i : 




0.000 0.005 0.010 0.015 

l«ool(eV) 

Figure 7: |m ee | distribution among matrices satisfying 2a angle and mass constraints 
described in the text, for N = 30 (blue, solid) and N = 100 (red, dashed). Also shown is 
the anarchical case (yellow, dot-dashed) under the same constraints. 

Finally, we have done a similar analysis for the distributions of CP phases; we have 
found that their distributions are relatively uniform, not providing us with any particular 
prediction. 

Other distributions 

To get a sense of how much our results depend on our particular choice of probability 
distribution, in Table [2] we look at the frequency of obtaining two large angles for a 
few different scenarios, which we take for the purpose of illustration rather than any 
particular physical motivation. We again require one angle with sin 2 (29) > 0.98 and 
another with sin 2 (26*) > 0.84. Table [2] also lists the percentage of N = 100 matrices 
which satisfy 0.004 < sin ($13) < 0.028 (2<j bound), after having first satisfied the 2u 
bounds on sin 2 (#12), sin 2 (#23), and the mass splitting ratio, described above. Case I is the 
distribution (5 = 1.1) we have already considered in detail, and is listed for reference. 

Case II demonstrates the importance of the second condition listed in Section 2- 
that enough right-handed neutrinos must be present so that each left-handed neutrino 



17 





5 


N = 3 


TV = 30 


N = 100 


N = 200 


Pioo(M 


Comments 


I 


1.1 


2.5 


10.1 


17.0 


19.0 


6 




II 


1.3 


2.5 


3.2 


6.7 


10.8 


11 




III 


0.9 


5.3 


17.6 


20.2 


20.9 


5 




IV 


0.9 


3.7 


16.3 


19.7 


20.5 


5 


No minimum cutoff 


V 


1.1 


3.1 


6.8 


15.4 


18.4 


7 


Hierarchical Majorana 



Table 2: Percent of events with sin 2 (2#)i argcst > 0.98 and sin 2 (2#) ncxt i argcst > 0.84 for 
different distributions and numbers of right-handed neutrinos, as described in the text. 
Also shown (denoted -Pioo(#i3)): for matrices satisfying 2a constraints for 9 12 , 9 2 3, and 
the squared mass difference ratio, the percentage which also satisfy the 2a bounds on 13 , 
with N = 100. 



is expected to have ^ 1 Yukawa couplings within a factor of a few of the upper bound. 
Here, the distribution steeply favors the low end of the spectrum (5 = 1.3), which makes 
sampling of the upper cutoff much rarer; as a result, many more right-handed neutrinos 
are needed to wash out the hierarchy. Whereas there is a clear enhancement by N = 30 
for 5 = 1.1, the rise in likelihood of large angles does not occur until much larger N in this 
case. Because the steeper distribution favors smaller mixing angles for a given number of 
right-handed neutrinos, this case performs somewhat better than the others for the 9± 3 
distribution. 

In case III, we see a distribution logarithmically skewed towards larger couplings (5 = 
0.9). This makes our large N effect much more dramatic, becoming noticeable for smaller 
values of N. Additionally, because this distribution is well behaved all the way to zero, 
it is not necessary to impose a lower bound on the Yukawa couplings for this value of 5, 
so we remove this requirement for case IV. Because only values close to the top of the 
distribution are important at large N, this only has the small impact of removing a bit of 
probability from the top. The change is more apparent at small values of N, where the 
lower parts of the distribution can be relevant in viable matrices. 

Finally, we turn to the choice of distribution for the Majorana mass matrix, which 
we have previously been taking to be anarchical and basis independent. In case V we 
take a very different approach, choosing the same hierarchical distribution (5 = 1.1) and 
cutoff structure for both Dirac and Majorana mass matrix elements. While this appears 
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to suppress washout effects somewhat, the impact is fairly weak, as the hierarchy is 
largely lost in the inversion of the Majorana matrix; only a small relative suppression of 
probability remains at N = 100. 

From these scenarios we see that while the quantitative behavior varies somewhat, the 
qualitative effect of increased probability of large mixing angles from many right-handed 
neutrinos applies broadly in the class of distributions we have laid out. Furthermore, 
we see that while there is some tension with a small value of sin(#i 3 ), all cases have a 
significant amount of probability lying in the experimentally favored region. 

4 Discussions and Conclusions 

In this work we have shown that the origin of the unique flavor properties of the neutrinos 
- small masses, and large mixing angles - may both have their origin in the seesaw 
mechanism. Even if the neutrino Yukawa couplings have large hierarchies, as are observed 
for the charged fermions, the hierarchical structure may generically become washed out if 
the seesaw involves a large number of right handed neutrinos, numbering perhaps in the 
tens or hundreds. We have given an explicit example showing the mechanism at work, 
with a statistical distribution of Yukawa couplings of a roughly scale invariant form, fit 
to the properties of the quarks and charged leptons. After integrating out the right- 
handed neutrinos, the probability of obtaining a mixing angle as large as the observed 
near-maximal value of #23, along with a second angle as large as the observed value of 
0i2 was found to be about 20%. The mechanism is fairly general, and may work with a 
variety of Yukawa structures, so long as various conditions are satisfied as described in the 
text. General predictions of the framework are a normal hierarchy for neutrino masses, 
and a neutrino mass matrix element |m ee | of about 1—6 meV. 

There are a variety of possible directions for future work. The first and most obvious 
one would be to try to build a top down model yielding hierarchical Yukawa couplings 
of an appropriate nature for our mechanism. This is non-trivial, since the most common 
approaches to producing hierarchical Yukawa couplings - flavor symmetries and extra 
dimensional wavefunctions - yield correlations across all of the couplings of a given field 
which do not lend themselves well to our framework. It may also be interesting to con- 
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sider the issue of leptogenesis j23j EI] in our scenarioj^j If the largest neutrino Yukawa 
couplings are of order one, as in the case of the quarks, then the seesaw scale we require is 
fairly large - of order 10 15 — 10 16 GeV. It then follows that washout processes will destroy 
any lepton asymmetry produced by the decays of the right-handed neutrinos. As a result, 
it may be necessary to move away from the thermal leptogenesis paradigm, and construct 
a model involving, for example, out of equilibrium inflaton decays to right-handed neu- 
trinos at temperatures much below the seesaw scale [26, 24J. One other possible future 
direction would be to consider the impact of anthropic selection effects on the types of 
Yukawa coupling distributions that we have been working with. One amusing possibility, 
somewhat orthogonal to the direction we have been pursuing here, is that the fundamen- 
tal Yukawa distributions for the charged fermions may actually be reasonably degenerate, 
but with strong anthropic selection effects constraining the first generation quarks and 
charged leptons [271 EE], and indeed, perhaps even the top quark (29], to be outliers. If this 
were the case, degenerate neutrino Yukawas might in fact be more representative of the 
fundamental Yukawa distributions than the hierarchical charged fermion ones. It might 
be interesting to see if this picture could be made to work at a quantitative level using 
some relatively degenerate ansatz for the Yukawa distributions, and putting anthropic 
constraints on various charged fermions masses. 
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